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We study the existence of monopole bound states saturating the BPS bound in N=2 
supersymmetric Yang-Mills theories. We describe how the existence of such bound states 
relates to the topology of index bundles over the moduli space of BPS solutions. Using an 
L 2 index theorem, we prove the existence of certain BPS states predicted by Seiberg and 
Witten based on their study of the vacuum structure of N=2 Yang-Mills theories. 
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1. Introduction 

Strong-weak coupling duality, or S-duality, has been conjectured in certain supersym- 
metric field theories and string theories. The existence of such a symmetry was originally 
proposed by Montonen and Olive for Yang-Mills theories P]. Among the testable predic- 
tions of S-duality is the existence of certain Bogomol'nyi-Prasad-Sommerfeld (BPS) bound 
states. For N=4 Yang-Mills, Sen verified the existence of such dyon bound states with mag- 
netic charge two . This provided a strong dynamical test for the existence of S-duality in 
this theory. For N=2 Yang-Mills theories, Seiberg and Witten have proposed a generalized 
S-duality involving the dependence of the theory on the Higgs field expectation value 0. 
For the case of N=2 Yang-Mills coupled to matter multiplets, they conjecture the existence 
of certain BPS bound states required by the singularity structure of the vacuum manifold. 
We propose to verify at least some of their predictions. 

Our approach to this problem involves quantizing the low-energy dynamics of N=2 
Yang-Mills coupled to matter. The existence of monopole bound states then reduces to 
a study of the spectrum of the Hamiltonian governing the low-energy dynamics. The 
existence of a bound state saturating the BPS bound is then equivalent to the existence of 
a zero mode for a twisted Dirac operator on the monopole moduli space. Since the moduli 
space is non-compact, we employ an L 2 index theorem to count the number of zero modes, 
and hence BPS states at magnetic charge two. 

We briefly summarize our results. N=2 Yang-Mills with Nf hypermultiplets has no 
BPS states at magnetic charge two for Nf < 3. For Nf = 3, we find two bound states: one 
with allowed electric charges 4n + 1 and the other with allowed charges 4n + 3, where n 
is an integer .0 These states are singlets of the SO (6) flavor symmetry. For the candidate 
S-dual theory with Nf = 4, we find bound states corresponding to the SL(2, Z) partners 
of the fundamental electrons; however, we also find partners for the heavy gauge bosons, 
implying that they exist as discrete states at threshold in this theory. 

In the following section, we derive the Lagrangian governing the low-energy dynamics 
of monopoles and dyons in super symmetric Yang-Mills coupled to matter. Section three 
describes the moduli space of BPS solutions, and the bundles of interest to us. The 
index computations needed to obtain the BPS spectrum are then presented. Section four 
provides a comparision with the states predicted by Seiberg and Witten. The final section 
is a summary and discussion. 

1 The electric charge of an electron is normalized to one as in S. 
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2. Supersymmetric Yang-Mills and the Collective Coordinate Expansion 

2.1. Pure N =2 Yang-Mills 

Collective coordinate expansions for N=2 and N=4 Yang-Mills around BPS monopole 
configurations have been described in |4|] . We require the slightly more general case of N=2 
coupled to general matter. Let us denote the SU (2) symmetry rotating the supersymmetry 
generators by SU (2)j. We shall proceed by first discussing the quantization of zero modes 
for pure N=2 Yang-Mills and then coupling to matter. In N=l superspace, the Lagrangian 
for N=2 Yang-Mills takes the form 

Lym = ^ J d*6&e 2V $ + j d 2 9TtW a W a + c.c^j , (2.1) 

where W a is a vector multiplet, and $ a chiral multiplet in the adjoint representation 
of the gauge group. In terms of component fields, the Lagrangian contains a gauge field 
A M , Higgs field cj), and an SU(2)i doublet of complex Weyl fermions rf . Note that gauge 
indices are suppressed for most of this discussion. Let us take the gauge group to be SU (2) 
which restricts any matter to the fundamental or adjoint representations. 

In the Coulomb phase, the gauge group is broken from SU{2) to U(l). The flat 
directions for the potential correspond to \(f), (f)'] = 0. Let <f> have vacuum expectation value 
(Tr (f> 2 ^ = ^f 2 which we choose to be large so that a semi-classical analysis is applicable. 
In the Coulomb phase, the theory possesses fundamental particles and solitons with masses 
saturating the BPS bound. The BPS spectrum consists of dyons, W-bosons, and with the 
inclusion of matter, electrons. Some of these particles may exist at the quantum level. As 
explained in 0, the mass of such a state is determined by the central extension of the 
supersymmetry algebra, and is given by 

M = V2\n e v + n m — ^— | (2.2) 
9 

in the semi-classical limit. Here n e ,n m are the electric and magnetic charges respectively. 
This formula remains true when matter is coupled since we shall assume all the electrons 
have vanishing bare masses. 

Since our interest is in checking for the existence of BPS states, we will eventually 
consider only low-energy fluctuations around a BPS solution. Such fluctuations are tan- 
gential to the moduli space of BPS solutions with charge fc, and so our computations will 
reduce to non-relativistic quantum mechanics on the moduli space. 
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2.2. BPS Field Configurations 

We shall follow the conventions of Wess and Bagger || . The supersymmetry trans- 
formations for the action (2.1) take the form: 

<L4 M = -ifjjGi^t 3 + iejfffjjf 

6<f> = V2e J T]j (2.3) 
8rf = a^F^e J - iV^a^D^ 

where we have set [(f), </>'] to zero, and where the e J a are anti-commuting parameters. We 
can further choose (ft to be real and search for static field configurations which preserve half 
the supersymmetries. Such configurations satisfy the BPS bound. In the gauge A Q = 0, 
with the magnetic field Bi = heijkF^ , we find the first order BPS equations 

B t = ±V2D i( ft. (2.4) 
which also follow immediately from the classical energy for a monopole configuration 




The second term is a topological invariant proportional to the magnetic charge and satu- 
rating the BPS bound (2.2). For magnetic charge k, a field configuration satisfying (2.4) 
depends on 4k parameters ip a . These 4k moduli are coordinates for an interesting non- 
compact hyperkahler manifold that will be discussed further in the following section. 
Associated to each modulus, <^> a , are zero modes (A l oa ,(ft oa ) in the expansion around the 
monopole configuration (A 1 , (ft). Such zero modes must be orthogonal to the directions 
generated by infinitesimal gauge transformations with compact support. This condition 
determines the adjoint- valued parameters ui a defined by: 

A^ a = d a A l + D'u a 

(2.6) 

<ftoa = {9a +U a )(ft. 

The standard technique used to quantize these moduli is to permit the f a to vary with time 
while omitting the zero modes from a perturbative expansion. We will further truncate 
the mode expansion of (A 1 , (ft) to the moduli dependent classical solution. 

We will expand the effective action to second order in the number of time derivatives 
and later to fourth order in the number of fermion fields. To ensure that the equations of 
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motion are satisfied to first order, the classical configurations must be corrected by terms 
of order one; specifically, the constraint A Q = must be modified to satisfy the gauge-field 
equations of motion where now 



D = <p a (d a + U a ) . 



(2.7) 



Taking into account such modifications, the bosonic piece of the effective action then 
describes a sigma model on the moduli space 



A more detailed discussion is given in [|J. 

2.3. Fermionic Zero Modes 

In the topologically non-trivial monopole background, the fermions possess zero 
modes. The inclusion of these zero modes in the effective action is readily accomplished 
by noting that the low-energy theory must be super symmetric since the BPS configuration 
preserves a single supersymmetry. The Callias index theorem states that the Dirac oper- 
ator for fermions in the adjoint representation in a monopole background of charge k has 
2k zero modes ||. We must therefore introduce 4k real Grassmann collective coordinates 
7 a in a mode expansion for the two Majorana fermions in the N=2 gauge multiplet. As 
in the case of the bosonic moduli, we allow the fermionic moduli to vary with time. By 
simply counting degrees of freedom, and using the pairing of bosonic and fermionic degrees 
of freedom required by (2.3), we can easily determine the effective action. The action is 
most simply expressed in terms of the superfields $ a = <p a + #7" where 9 is a Grassmann 
coordinate |7j: 



The super-covariant derivative, D, is given by D = —i-§g+8-§i- This sigma model possesses 
more supersymmetries than are apparent from this superspace formalism since the mod- 
uli space is hyperkahler, but we will not need to make those additional transformations 
explicit. Quantization of such a quantum mechanical model is a well-studied problem. 
Supersymmetric ground states correspond to zero modes of the Dirac operator on the 
moduli space, or equivalently to anti-holomorphic closed forms on Mk since the space is 
Calabi-Yau || Q . We can now move easily to the case with matter multiplets. 





where 




(2.9) 




(2.10) 
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24- N=2 Yang- Mills Coupled to Matter 

The Lagrangian for an N=2 hypermultiplet contains two chiral N=l superfields M 
and M in conjugate representations of the gauge group. With bare masses set to zero, the 
matter Lagrangian is given by: 



L M = J d A 9 (M ] e 2V M + M^e~ 2V M^j + (^2 J d 2 6M$M + c.c^j . (2.11) 

The components of the hypermultiplet are an SU(2)j doublet of scalar fields {m^rfv ) to- 
gether with complex Weyl fermions ^A, A j or equivalently two Majorana fermions. We 
shall discuss the case where M is in the fundamental representation of SU(2) with her- 
mitian generators T . In the Coulomb phase, the potential energy has no additional flat 
directions. Further, there are no zero modes for the scalars in the monopole background 
since the operator 

-D i D i + (\/2(j) l T 1 ) 2 , 

is positive. However, the fermions do possess zero modes which form a bundle over the 
moduli space M^. The bundle has dimension k and transition functions in O(k). For each 
Majorana fermion zero mode A on , we introduce a Grassmann collective coordinate t(j n , 
n = 1, As usual, the ip n are dynamical and vary with time. The fermionic kinetic 

term for each Majorana fermions in the Lagrangian (2.11) provides a natural bundle metric 
hnm(^p) = / d 3 x \on(x, Lp)"f° \o m (x , (p) for the low-energy theory. The other ingredient 
needed to fully describe the effective action is a connection for the bundle. The connection 
also follows from the full action where we recall that Aq can now contain terms with two 
fermions as well as the terms of order one in time derivatives. Let Dq be the restriction 
of D obtained by setting the fermions (7, ip) to zero. An expression for the connection O 
then follows from the term 

/ d 3 x Aon7°A)Aom = / A on7 V (da + u l a T l ) X 0m 

where u a satisfy (2.6). The effective action can then be described using fermionic super- 
fields = ip n + OF 71 where the second component is an auxiliary field, and 

S eff = J dtd9g ab {^D^ b + ih nm {^ na {D^ + D^ a ^l). (2.13) 



The additional index a affixed to the fermionic superfields runs from 1, . . . ,2Nf since 
there are 2Nf Majorana fermions in the fundamental of an SO(2Nf) flavor symmetry 
described in . The symmetry properties of the bound states under this flavor group will 
be discussed further in section 4. The supersymmetry generator for this theory corresponds 
to the Dirac operator coupled to the 0(k) bundle connection. Supersymmetric ground 
states then correspond to normalizable zero modes of this twisted Dirac operator. 

The case that we will focus on for the remainder of the paper is magnetic charge 
k = 2. Since the connection has an abelian component in this case, states in the Hilbert 
space of this theory are labelled by their U(l) charge. Charge conjugation symmetry pairs 
a state with £7(1) charge n to one with charge -n. In the case of Nf hypermultiplets, the 
spectrum of U(l) charges takes the form 

\n>,\n+l >,..., \n + 2N f > (2.14) 

where n is the U(l) charge of the Fock vacuum. Charge conjugation then implies that 
n = —Nf. To summarize: the task of finding BPS bound states at magnetic charge two is 
equivalent to that of finding normalizable zero modes of the twisted Dirac operator acting 
on states with U (1) charge —Nf, . . . , Nf. 



3. An Index Theorem for the Atiyah-Hitchin Manifold 

To count the number of normalizable zero modes of the twisted Dirac operator, we 
will employ index theory. Clearly, we need to understand the two-monopole moduli space, 
its metric, the bundle and its connection. Fortunately, these topics have been well-studied, 
and most of the structures we require are known. We will benefit greatly from the work 



of ||10|| , ||11|| and |L2| in this analysis. However, the index computation is subtle since the 
moduli space is non-compact. To count the number of L 2 modes will require a careful 
analysis of boundary effects. First however, we must describe the geometry of the moduli 
space. 

3.1. Monopole Moduli Space 

As discussed in section 2.1, the BPS equations (2.4) for magnetic charge k have families 
of solutions with 4k parameters. For k = 1, these parameters describe translations and 
"large gauge transformations" of the standard BPS monopole solution. The moduli space 
in this case is just R 3 x S . For general k the space has the form 



6 



Mk = R 3 x (3.1) 

where M® is a 4k — 4 real dimensional hyperkahler manifold equipped with a hyperkahler 
metric [ 10 1 . It is often useful to consider the k— fold cover M k = R 3 x S 1 x M%. The metric 



on Mfc is flat in the R 3 x S 1 factors, so we will focus our discussion on M% - specifically, 
Ml 

Invariance of the BPS equations under the Euclidean group amounts to, in part, an 
SO (3) group of isometries on M®, the generic orbit of which is three dimensional. We 
identify so(3) with the left invariant vector fields on SO (3) in the usual manner. Let 
{Xi, X2, X^} be a basis for so(3) satisfying 

[X ll X 2 \ = —X 3l [X 3 ,Xi] = — X 2l [X 2 ,X 3 ] = — Xi, 

and let {0^} denote the dual frame. They satisfy the relation doi = \tijk&j A ak- 
in this notation, the metric on is constrained to be of the form 

ds 2 = f(r) 2 dr 2 + ai (r) 2 a 2 + a 2 (r) 2 a 2 + a 3 (r) 2 a 2 . (3.2) 

Where convenient, the Oi will be described by 50(3) Euler angles < 6 < tv,0 < <fi < 
2tt,0 < V < 2tt : 



o"i = — sin tpd9 + cos ijj sin 
o"2 = cos tftdO + sin ip sin 9d<p 
03 = dip + cos 6d(j). 



also has the identification 

(r,e,M) = {r,ir-e,<l> + Tr 7 -il>). (3.3) 

Anti-self-duality of the curvature - following from hyperkahlerity - tells us that 

2a2&3 da± , , 2 2 

as well as the equations obtained by cyclically permuting a\, 0,2, 03- 

By redefining the radial coordinate if necessary, we can take / = — az/r, and the range 
of r is then from n to 00. For this choice the large r dependence of the metric is found to 

be 

/T~2 

ai ~ a2 ~ rW 1 



up to terms suppressed by e r NTT]]. Near the bolt coordinate singularity at r = %, we have 



a i 


« 2(r 


-7T) 




a 2 


1, 
~ 2 V 


-7T) 


+ 7T 


a 3 


1, 
~ 2 V 


-7T) 


— 7T 



As discussed in the previous section, the zero energy solutions to the Dirac equation 
form an O(k) bundle over the moduli space The bundle is trivial over the R 3 factor in 
the decomposition of the moduli space given by (3.1). Further, Hitchin has shown that the 
connection on this bundle has anti-self-dual curvature.i Let us first consider the case k = 1 
where the index bundle is an 0(1) bundle. While obviously flat, the bundle is not trivial, 
and is once twisted over the S 1 factor. We now focus our attention on the nontrivial case 
k = 2. Let Ind2 denote the restriction of the bundle to (S 1 x M®)/^- The obstruction to 
orienting the bundle arises from the explicit discrete Z2 identification in (3.1), which we 
call Is- In coordinates, Is is the identification 



(X; r, 9, </>, = (x + tt; r, 9,(f),^ + 7r), (3.4) 

where < x < 27r is a coordinate for S . We can therefore pull back this bundle to 
S 1 x M® and make a choice of orientation giving a (7(1) bundle Ind2- 

The anti-self-duality of the bundle curvature, and its invariance under the SO (3) 
subgroup of the Euclidean group of symmetries of the BPS equations, fix the curvature to 
lie entirely in and have the form 



O = da A o"i + ao"2 A cr 3 . (3.5) 
The function a(r) falls off as e~ r l 2 as r — > 00, and has a normalization a(n) = ±i [12 



The sign ambiguity depends on the choice of orientation, and will play no role in the 
following analysis. 



As referred to in [12]. 



3.2. The Index Formula 

In this section, we consider the restriction of Ind2 to (p, M®) = M® , where p is any 
point on S 1 ^ I inherits a connection from Ind2- We wish to compute the dimension of 
the kernel of the Dirac operator D on spinors with values in I n . First we need a vanishing 
theorem. 

Proposition 3.2.1 Let M be an infinite volume spin four-manifold with zero scalar cur- 
vature. Let E be a line bundle over M whose curvature r E is an anti- self- dual two-form. 
Then the kernel of the Dirac operator D% on S~ <8> E is zero. 

Proof. Integrating by parts and using the Bochner-Lichnerowicz formula ( ||13|| , Theorem 
8.17) we have 

ll^/H 2 = l|V/|| 2 + ( J R E /,/), 

where R E = J2i<j e i e j r F,ji an d e i denotes Clifford multiplcation by the i th vector in a 
frame. In an oriented frame, we may write the projection onto S~ ®B as (1 — eoeie2e3)/2. 
The anti-self-duality of r E implies 

} j e i e j r E j (l - e e 1 e 2 e 3 )/2 = 0. 

i<j 

Hence, 

ll^/ll 2 = liv/|| 2 , 

for / a section of S~ ®E. An element of the kernel of is therefore covariantly constant. 
This implies its norm is covariantly constant. On an infinite volume manifold the only 
covariantly constant square integrable section is zero. □ 

Corollary 3.2.2 The L 2 index of the Dirac operator 

D+ : L 2 (M 2 °, S+ ® I n ) -> L 2 (M 2 °, S~ ® I n ) 

is the dimension of the kernel of D + . 

Thus we are left with the computation of the L 2 index. We recall the modifications 
required to compute the index on a complete noncompact manifold. Let e~ sD denote the 
heat operator associated to D 2 . As s — > oo, e~ sD converges weakly to projection, IT, onto 

3 Though Ind2 wraps nontrivially around the S 1 , there is no ambiguity in defining /, up to 
isomorphism. 
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the kernel of D. Hence, the integral of tre~ sD (x, x) over a compact subset C converges 
to J c trU(x,x)dx, and J c trre~ sD (x,x)dx — » J c trrU(x,x)dx, where r denotes Clifford 
multiplication by the volume element, r = eoeie2e3 in the notation of the preceding 
proposition. The L 2 index of D is given by the integral 

Ind(-D) = / trrll(x, x)dx = / Zzm s ^ 00 trre _s ' D (x,x)dx. 
Jm° Jm° 

On a compact manifold, Y, with Dirac operator Dy, J Y trTe~ sDY (x,x)dx is independent 
of s; hence, in the above discussion the s — > oo can be replaced by s — > 0, and one obtains 



Ind(-Dy) = lim s ^Q J trre sDY {x,x)dx. 



In the latter limit, the traces are easily computed. In general -^trre sD (x, x) can be 
expressed as the divergence of some vector field V(s), with 
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divV(s)dx = J tre^rDe sD da. 

C JdC 

Here eo is Clifford multiplication by the unit outward normal to dC and da is the induced 
volume form on dC. (See ||.) This explains the s independence of the trace in the compact 
case and yields the following expression for the index on (see, e.g. [0, |jl5| , and 0): 



Ind(-D) = limL^ooilirris^o I trre sD dx + / tre^rDe sD dads). 

Jr<L JO Jr=L 

The Atiyah- Singer index theorem gives 

limL^ooilirris^a / trre~ sD dx) = 

Jr<L 



24 



-i-j / Tr(RAR) + -L[ Tr(fMfi). 
• 8tt 2 J m o 8tH J m o 



We write 



Ind(D) = — [ Tr(R A R) + -L I Tr(OAfi)+5 D , (3.6) 
24 • 8tt z J m o Htt z J m o 



where 



5d = HrriL^oc / / treorDe sD dads. 

Jq Jr=L 
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From the form of the metric (3.2), one computes the curvature in the orthonormal 
frame 

{v , v\,V2, ^3} = {fdr, 010-1,(120-2,030-3} 



to be 



-Rio 


= R 


R20 


= R: 


R30 


= R 



23 = ^jO:A)u Au 1 + ^- 3 (A + B + C-l-2BC)u 2 A^ 
31 = ^j(i:B)u Au 2 + ^- 3 (A + B + C-l-2AC)u 3 Au u 
12 = i(|C>oA^3 + i(^ + fi + C-l-2i%A, 2) 



where 



A = (a 2 + a l ~ a 2 )/2a 2 a 3 , 
B = (a 2 + a\ - a|)/2a 3 ai, 
C = (a 2 + a\ — a\)/2a\a 2 . 
Anti-self-duality holds for the choice of orientation u A v\ A v 2 A z/ 3 , and we see that 

Tr(R A R) = d [{-4(A + B + C - l) 2 + 16 ABC) 01 A 2 A 03] , 

which by Stokes' theorem (note that the boundary orientation induced here is — v\ Av 2 Az/3) 
and the asymptotic formulas for a±, a 2 , and a 3 gives 
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If 1 

^ / Tr(RAR) = . 

•8W M o 1 J 6 



Note that we have divided by a factor of two due to the Z 2 identification (3.3). 

The next term in the bulk contribution is also simple to compute, since we have, from 
(3.5), 

O A O = n 2 d(a 2 A 01 A 2 A 03) 
for the bundle I n . Using ct 2 (oo) = and a 2 (n) = 1/4, we find 

At / Tr(fiAfi) = — . 

Combining the terms, we arrive at 

Ind(D) = n 2 /8 - 1/6 + 8 D . (3.7) 
We now turn to the calculation of the boundary contribution So- 
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3.3. An Equivalent Index Problem 

For K some large positive constant, equip (K, oo) x 50(3) with the metric f(r) 2 dr 2 + 
J2 a l( r ) a i- This metric descends to a metric on (K, oo) x SO(3)/Zi2, where the Z 2 action 
is the one generated by (3.3). This is just the Atiyah-Hitchin metric described in (3.2) 
near 00. 

In order to reduce the computation of the index of the Dirac operator to computations 



essentially the same as those carried out in JI3J], we will first make a conformal change in 



the metric of the moduli space in a neighborhood of 00. This change is not essential to 
computing the defect 5d but simplifies the estimation of error terms involved in construct- 
ing approximations to resolvents and heat operators. By the conformal invariance of the 
Pontrjagin classes, we know that such a conformal change will not change the value of the 
bulk contribution. Hence, in order to justify it, we need only check that it preserves the 
index. First let us specify a new metric which, for large r, is given by 



fir? , 2 . a 2 



-dr 2 + J2 



* 

2 a i 



2 1 2 (3.8) 

dv 2 4 

= ~Z + (1 + A M + °l + ( r 2 (1 _ 2/r) 2 + 

with Aj doubly exponentially decreasing functions of t (i.e. 0(e kt e~ e )). Set t = ln(r). 
The metric is then 



dt 2 + (1 + A x )a 2 + a 2 + ( e2t(1 _ 2e - t)2 + As)a 2 . 

We fix an orthonormal frame {Yi} on the Z 2 cover of our space with Yq = J^, Yi = 
(1 + A 4 )X 1 , Y 2 = X 2 , and Y 3 = (e t /2 - 1 + A 5 )X 3 , with A 4 , A 5 doubly exponentially 
decreasing. Of course this frame does not descend to one defined globally on the Z 2 
quotient, but that will not affect our computations. Quantities involving the squares of 
these vector fields will descend to the quotient. 

Neglecting doubly exponentially decreasing terms, we have the commutation relations: 

ln,Y 2 ]~- T ^Y 3 , 

[y 3 ,n]~-e*(i/2-e- t )y 2 , 
[y2,y 3 ]~-e*(i/2-e-*)n, 

[F ,y 3 ]~(l-2e- t )- 1 F3. 
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From these relations, we can check modulo doubly exponentially decreasing terms that 

(Vy.Fi, Y 2 ) = -(e*/2 - 1 - (e*/2 - l)" 1 ^), 

(Vy 1 y 2 ,y 3 ) = -(e*-2)- 1 , 

(Vy 2 y 3 ,F 1 ) = (e t -2)- 1 , 
(Vy 3 y ,l3) = -(l-2e- i )- 1 . 

For the remainder of this section, let D denote the Dirac operator with respect to 
this new metric and denote by D the Dirac operator associated to the old metric g = 
e 2ioga 2 g_ reca u ( see J|3j ) that there is a local isometry between the spaces of spinors 
determined by the two conformal structures such that 

Ds = e- 3lo9a2/2 D(e 3lo9a2/2 (p(s)). 

Hence, the map 

$ : 3 -> e 3 ^° sa2/2 (/)(s) 

takes harmonic spinors to harmonic spinors. This would induce an isomorphism between 
spaces of L 2 harmonic spinors if $ preserved the L 2 condition. 

Proposition 3.3.1 The map $ is an isomorphism between spaces of L 2 harmonic spinors. 
Proof. Let h G L 2 (g). Then 

oo > J \h\ 2 dv- g ~ J \h\ 2 e 4t dv g . 

Hence, J \e 3t / 2 (j)(h)\ 2 e t dvg < oo, and $ maps 1? harmonic spinors to 1? harmonic spinors. 
To obtain an isomorphism, we must also show that the inverse map also preserves the L 2 
condition. It suffices to show that for H G Ker(Z)), e i ^ 1 / 2+a ) H G L 2 (g) for some positive a. 
We shall prove this estimate below. First we need a preliminary discussion of the holonomy 
of the index bundle, I n . 

The fundamental group of SO(2>)/7i2 is generated by the inclusion of the circle K 
obtained by exponentiating the Y3 vector. This is the fiber of the fibration 

{L} x 6X>(3)/Z 2 -> {L} x SO(3)/KZ 2 . (3.9) 

According to |T^], the holonomy about this circle fiber is given by e ±l7rri / 2 when L = 00. 



They show this by integrating the curvature over the Atiyah-Hitchin cone [1C]. Hence, for 
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finite L, the holonomy differs from this factor by e l7re ( L ), where e is doubly exponentially 
decreasing. This can be seen by integrating the doubly exponentially decreasing curvature 
over the subset of the cone with r > L. Note that the Z 2 action prevents us from globally 
fixing the sign since Y3 is only globally defined up to a factor of ±1. We will assume Y3 
chosen so that the holonomy ~ e -™n/2+ie{L) _ 

It is more convenient to work with periodic frames than with covariantly constant ones. 
Multiplying a 13-covariantly constant frame for I n by e _r ^( n / 4_e / 2 ) gives a connection form 
of ^Y" 3 * modulo a doubly exponentially decreasing term. 

Taking now a frame for the spinors determined by our (periodic) frame {Yq, Y\, Y2, Y3} 
and the periodic frame for J n , we can Fourier expand sections in the Y3 direction. This 
can even be done globally, although the Fourier coefficients then become sections of a Z 2 
quotient of powers of the Hopf bundle over S 2 . 

Setting T = Y , we can write 

D 2 = —T 2 - V 2 n - V 2 Y2 - V 2 Ya + Vv y ,y„ 

modulo rapidly decreasing terms. Then taking the k th Fourier component with respect to 
the circle action associated to Y3, we have 

D 2 = — T 2 + T — (e*/2 - l) 2 (ik + in/4 - 1/2(1 - [e*/2 - l}- 2 /2) ei e 2 
~ 4(e*/2- l) 2 6063)2 

plus a positive operator and rapidly decreasing terms. This operator is conjugate to 

- T 2 + 1/4 - (e*/2 - l) 2 (ik + in/4 - 1/2(1 - [e*/2 - \\~ 2 /2)e x e 2 
e t 

~ 4(e*/2- l) 2 6063)2 

plus a positive operator and rapidly decreasing terms. The smallest eigenvalue of — (ik + 
in/4 — l/2eie2) 2 is 1/16 when n ^ 2 (mod 4), since k is an integer. Hence, a standard 
maximum principle argument (or differential inequality) says that for n ^ 2 (mod 4), any 
L 2 harmonic spinor is doubly exponentially decreasing and therefore maps to an L 2 (g) 
harmonic spinor under the previously described map. Moreover, this estimate implies that 
D is Fredholm with no continuous spectrum in the case of restricted n. We are left to prove 
the proposition in the case n = 2 (mod 4). In this case, (ik + in/4 — l/2eie 2 ) can have a 
kernel. On this kernel, the —V 2 term is dominated by ~( 4 ( e t / f 2 -i)^ e o e 3) 2 = 1/4 modulo 
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exponentially decreasing terms. Hence, the decay of the component corresponding to the 
kernel of (ik + in / A — l/2e\e2) is, by a maximum principle argument, of the order 0(e~ t ^ 2 ). 
As the volume of the SO (3) orbit at t is 0( e -*), e^+^H G L 2 (g) for H G Ker(D). This 
implies that <E> _1 takes L 2 harmonic spinors to L 2 harmonic spinors, completing the proof 
of our proposition. □ 

This proof shows that the essential spectrum of D 2 is bounded away from zero. Hence, 
D is Fredholm. A similar argument works to show when n ^ 2 (mod 4) that D is Fredholm. 

3.4- Parametrices 

For the convenience of the reader, we gather in this section aspects of the construction 
of the approximate heat kernels. The idea here is to determine an operator which is an 
appropriate approximation to e~ sD for use in computing 5d- We will define such an 
operator and find that only a few relevant components contribute to our calculation. 

Using the functional calculus, we write 

e~ sX (D 2 — X)~ 1 d\, 

where 7 C C is a simple curve surrounding the spectrum of D 2 . Thus, it suffices to 
approximate (D 2 — A) -1 . 

We construct the semilocal approximation to (D 2 — X)~ 1 inductively using a continuous 
Fourier transform in the base variables (in a coordinate neighborhood) and a discrete 
Fourier expansion in the fibers. More precisely, we consider contractible open subsets of 
RP 2 (the base of the fibration) over which the fibration is trivial. As in section 3.3, we 
Fourier expand the sections of the twisted spinor bundles. Fix a coordinate neighborhood 
V on the base space of the fibration. On the k th Fourier component D 2 has the form 

-T 2 +T-X7 2 Yi - V 2 Y2 - (e*/2 - l) 2 (ik + in/ A - 1/2(1 - [e*/2 - l]" 2 /2)e l e 2 

~ 4(eV2-l)2 e ° e3)2 

plus rapidly decreasing terms. 
Let 

K 2 = (e*/2 - l) 2 (ik + in/4 - 1/2(1 - [e*/2 - l]" 2 /2)e 1 e 2 ) 2 '. 

Let || v || 2 = ^2g^Vij. In a nice frame, let / = J2 k fk{x)e lk ^ denote the Fourier expansion of 
a section /, and let fk(v) denote the Fourier transform of Then the action of D 2 on / is 
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given via the inverse Fourier transform by operating on //. by multiplication by \\v || 2 + K 2 
plus some lower order operator L, which is readily computed. Hence, to invert (D 2 — A) 
approximately, it suffices to invert (\\v \\ 2 + K 2 — A + L) approximately. We construct such 
an inverse in the form Ylf=o (\\2-Kv\\ 2 +k 2 -\y+ 1 • ^ or ^ ne S ener i c case, K is very large and 
so high powers of (||27™|| 2 + K 2 — A) -1 are rapidly decreasing. Hence, the numerators are 
constructed by an inductive process, so that when acted on by D 2 — A (in the guise of 
|| v || 2 + K 2 + L) we obtain 1 plus a high power of (||27rv|| 2 + K 2 — A) -1 . The inductive 
construction is given as follows. Set 

ho(x, x', A, v, k) = Identity. 

Write 

(D 2 - X^n^e^-^e^-^ipnvW 2 + K 2 — A)"'" 1 /^, x', A, v, k) 
= {2Ti)- 1 e lk ^-^e l2 ^ x - x '> v {\\2Tiv\\ 2 + K 2 - X)- l hi(x,x', A, v, k) — 
2V((2Tr)- 1 e ik ^-^e i2 ^ x - x '> v ) • V{(||2t™|| 2 + K 2 — Xy^h^x, x' , A, v, k)} 
+ (2n)- 1 e lk -^-^e l27T{x - x ' ) - v A 2 {(\\2nv\\ 2 + K 2 — X)- l ~ l hi{x, x', A, v, k)} 
= (2n)- 1 e lk(i '- i '' ) e l27T ^ x - x ' ) - v (\\2nv\\ 2 + K 2 - X)- l [hi(x,x', A, v, k) + Ri], 

where A 2 = D 2 - K 2 . Set 

h l+1 = -(R^e^-^e^-^ipTvvf + K 2 - X) l+1 . 



Then one obtains formally 

N 



n2_uVV ! (n^-i^k-g,-^') r27,(x-x')-v hi(x,x',X,v,k)dv 
rhJ*.* (\\2nv\\ 2 + K 2 - XY+i 



Identity + / (Rn)cIv, 

and we take for some large N, 



N 



^to'^ (\\2nv\\ 2 + K 2 - X)^ 

for our semilocal approximation to (D 2 — A) -1 (pre- and post-multiplied by cutoff functions 
in the usual way and summed over a cover, etc.). 
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Write 

Dhi(x, x, A, v, k) = ^(||27rt;|| 2 + K 2 - X)- a hi, a (x, v, k), 

a 

and 

h u<J {x,v,k) = ^,^,B,c(^)f A (e*/c) B e- tC , 

0<\A\+B<2a+l,C>0 

with hi yay A,B,c(x) bounded. 

In the following section, we shall study traces of the operator 



o J-i 



_ sX . ^, (T ,A,B,c(^)^ A (e i fe) g e tc )dvd\ds 

' R 3 (\\2>irv\\ 2 + k 2 - xy+i+v 



(( ke t)[B+\A\-2l-2*-C+l]) < 0((/ ce *)[2-2/-C]_ 

Hence, we see without yet using the trace identities that only the terms with / = or 1 
can contribute to the trace. Let us now recall the basic trace lemma. 

Lemma 3.4.1 For • • ,ij distinct and j > 0, 

tre h ■ ■ ■ e ij = 0. 

So in computing the trace of eorDhi(x, x, A, v, k) we get a nonzero contribution only from 
those terms which introduce enough Clifford factors to cancel e$r = —e\e<iez - and no 
more. No term in Dh contributes a factor of e3 without also contributing an extra eo 
except the term 

e 3 (e*/2 - l)(ik + in/4 - 1/2(1 - [e*/2 - l]" 2 /2)e 1 e 2 )h 

so this is the only term which can have nonzero trace. We will compute the contribution 
of Dho to this trace. A similar computation shows that the contribution of Dh\ is rapidly 
decreasing. 

We compute 

, v ho(x, x', A, v, k)dvd\ 
2m J" ^W R 3 V "' V " ' (||2t™|| 2 + X 2 - A) 



— i e~ tx Y I (2n)- 1 e lk ^' ) e l2 ^ x - x ' ) 



(il>-il>') p i2n(x-x')-v dvdX 



2m j 1 V^ R3 (||2t™|| 2 + K 2 -A) 

:J2(^t)- 3 / 2 e-\ x - x '\ 2 e- tK2 (2n)- 1 e^-^\ 
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Recall we need to compute the trace of 



J2 e ore 3 (e t /2 - l)(ik + in/4 - 1/2(1 - [e*/2 - l]" 2 /2)eie 2 )(47rt)- 3 / 2 



x 



A; 




along x = x', ijj = ip' where we recall that r is the volume form eoe^es. Integrating this 
trace over the S 1 fiber reduces us to computing 



£Veiea(e*/2 - + in/4 - 1/2(1 - [e*/2 - l]" 2 /2)eie 2 )(47rt)- 3/2 e - tK2 . (3.10) 



5.5. Defect Computations 

We can now turn to the computation of <5d defined in (3.7). Let denote Clifford 
multiplication by Yj,. In our situation, we have the following 1? index theorem as we have 
previously discussed. 

Proposition 3.5.1 The index on I n is given by the expression 



the special case of a smooth divisor. Hence we will restrict ourselves here to providing 
an outline of the methods involved. The computation requires two steps: first we need to 
construct an explicit approximation E s to e~ sD so that 



Jt=L Jo 

Using the results of the previous section, E s takes the form (suppressing patching and 
partitions of unity etc.) : 




The computation of So is essentially the same as the index computation in in 





x 



hi(x, x', A, v, k)dv 



(\\2irv\\ 2 + K 2 - XY+ 1 ' 



Second, we need to compute 
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For N large and n ^ 2 (mod 4), the presence of the product of — (e*/2 — l) 2 and the 
square of the Fourier coefficients in the exponent can be used to show that D(E S — e~ sD ) 
has very small trace norm when multiplied by a cutoff function supported near oo. 

When n = 2 (mod 4), this construction does not yield good error estimates on the 
kernel of (ik + in/ 4 — l/2eie2). On this subspace, we must use a parametrix construction 
which is global on the entire SO (3) orbit. This is similar to the construction of the 
parametrix for singular m in [13] - see section 6 in particular. There is a natural limiting 
Dirac operator as L — > oo on this subspace, and we construct our heat operator (restricted 
to this subspace) as a perturbation of the heat operator associated to this Dirac operator. 
We need not compute this operator explicitly. It suffices to note that it does not contribute 
any Clifford multiplication factors to prevent cancellation in the super trace. The necessary 
Clifford factors arise in the perturbation expansion, but with 0(e~ L ) coefficients. Hence 
this term contributes nothing to our trace. The trace over the kernel of (ik+in/4— l/2e\e-z) 
of the parametrix obtained by the semi-local construction also contributes zero. Hence, 
we may use the semi-local construction for our entire computation. 

Set U = (e L /2 — Returning to the computation of the traces with the semi-local 

parametrix, we see from (3.10) that 

f /"OO /"OO 

/ / tre rDe~ sD2 ds = (4tt) V / tre x e 2 V~ x {ik + in/4 - (1/2- 

Jt=L JO k JO 

t/ 2 /4)eie 2 )e^ (lfc+iri/4 - (1/2 - c/2/4)eie2)2 x 
(4ns)- s / 2 ds + 0(e- L ). 
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We simplify and apply the Poisson summation formula to write this as 

/*oo p 

An J2 tre 1 e 2 U- 1 (ix + in/4:- (1/2 - £/ 2 /4)eie 2 ) x 
Jo p Jr 

e ^( l x +l n/4-(l/2-C/ 2 /4)e 1 e 2 ) 2 (47rt) -3/2 e -2^x M + q^- L) 

= 4tt J2 tre 1 e 2 U- 1 ixe^ x2 (4nt)- 3/2 e- 2lTipx x 
Jo p Jr 

e 2^ P (n/4 +l (l/2-[/ 2 /4)e ie2 ) rfa ,^ + q^-L^ 
poo 

= 4ntr / Ve 2 ^/ 4+ ^ 1 / 2 - c/2 / 4 ) eie2 )e 1 e 2 ^- 1 (-27rz)- 1 x 
Jo p 

d_ n 2 p2u2/t ^ /2ut _ 1/2 3/2dt + o(e _ L) 
op 

= try e 2 ^/ 4+ ^ 1 / 2 -^/ 4 ) eie2 )e 1 e 2 z^- 1 (-27rz)- 1 |-^^- + 0(e- L ) 
^ " op 2-K 2 p z U 

= trY e 2^p(»/4+i(i/2-t/ a /4)e iea ) eie 1 +0 (e- L ). 

We now eliminate terms using the trace identities. Replacing the factor e ^p(u 2 /2e 1 e 2 ) j-,y 
e irp(u 2 /2eie 2 ) _ ^ ^ oes not c h an g e the trace, as treie 2 l = 0. Moreover, we have that 

e n P (u 2 /2e ie2 ) _ 1 = np {ij^ /2 ei e 2 ) + 0(e" 4L ). 

Hence, the defect reduces to 

tr^e^ ( _ 1) P^2 eie2eie2 _l +0 {e-% 

which finally gives, 

5 D = J2e ntpn/2 (-l) p+1 1 



7T 2 P 2 

p^O F 

We compute this for different values of n (mod 4). When n = 0, 

S D = C(2)/tt 2 , 

where ((s) denotes the Riemann zeta function. As is well known, £(2) = 7r 2 /6. Hence, for 
n = (mod 4) 

5 D = 1/6. 
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If n = 1, (mod 4) we obtain 

-TT" 2 ^p-2 e - P 3/2 = _ n -2 J2p- 2 (-l)PtP = 1/2 ^7T-V 2 (-1) P+1 = 1/24. 

p/0 p/0 p>0 



If n = 2, (mod 4) we obtain 



-7r- 2 ^p- 2 e^ 2 = -l/3. 



If n = 3, (mod 4) we obtain 



-TT" 2 



^-2^5/2 = ^-a^ya-p = _i/ 2 ^ 7r - 2 p- 2 (-l) p = 1/24. 

p/0 p/0 p>0 




Before proceeding, let us summarize the results of this index computation: 



dimKer (£>/«) = 7178 - 1/6 + 



3.6. Sections o/Ind2 and Determining the Electric Charges 

Until now we have been working on the space M 2 °, although we are actually interested 
in sections of the index bundle Ind2 over (S 1 x M 2 )/I^, where ^3 is a Z 2 involution. Since 
the bundle is trivial over the R 3 portion of the monopole moduli space, we can restrict 
our attention to the dependence of the bound state wavefunction on S 1 x M 2 . To obtain 
the full wavefunctions, we can simply multiply by a factor of e ip ' x , giving the particles 
momentum. 

To define sections of the quotient bundle Ind2, one must equivariantly lift the Z 2 
action from the space S 1 x to the total space of the bundle Ind2 over S 1 x M^. The 
M® piece of this bundle has been discussed in section 3.1. Since Mi? is homotopic to S 2 , 
it is convenient to describe this piece of the bundle as a bundle over S 2 - it is in fact the 
Hopf bundle (the total bundle over M 2 ° is then isomorphic to the pullback) |12[ . The Hopf 
bundle S 3 — > S 2 is the quotient map z = (zi,z 2 ) 1— > \z\,z 2 \, where \z\\ 2 + \z 2 \ 2 = 1 and 
[zi, z 2 ] is the point in CP 1 = S 2 denoting the complex line through (z\, z 2 ). The fiber is 
clearly £7(1). Then by the homotopy, we can think of the total space of the bundle Ind2 
over S 1 x M 2 as a bundle P over S 1 x S 2 , defined as follows |fl2fl : 



P = Rx S 3 / {(t, z) ~ (t + 27r, —z)} . (3.11) 
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The bundle map is trivial on the first coordinate and the Hopf map given above on z. 
Note that on the base t ~ t + 2n, but we don't have (i, z) ~ (t + 2n, z) on the total space. 
Note, too, the nontrivial S 1 twist encoding the holonomy over the S 1 factor. The bundle 
over M2 given by pulling back the Hopf bundle in the manner discussed above is just the 
bundle / of the previous sections. Now the relation (3.11) is equivalently 

rnd 2 = Rx //{(*, e) ~ (t + 2tt, -e)}, (3.12) 

where — e is — 1 ■ e, with —1 acting along the fiber. 

To get to Ind2, we must quotient this bundle by Z 2 . This procedure involves lifting I 3 
to Is acting on the total space of Ind2- By "lift," we mean a commutative diagram 

h 

P — P 



h 

S 1 x S 2 — > S 1 x S 2 . 

In describing the lift, we will for convenience work with the bundle P defined in (3.11) and 
use the same notation I3 and ^3 - our statements can be "pulled back" to Ind2- In terms 
of S 2 , the -Z3 on acts downstairs as the antipodal map on 5" 2 : (2:1,2:2) | — ► (—Z2,zi)- 
Then J 3 maps 

(t, (z 1 ,z 2 )) -> (t + ir,(-z 2 ,zi)). 

One easily checks that ^3 squares to the identity as a result of the equivalence (3.11), and 
that the diagram is commutative. Now the quotient 

Ind 2 = Ind 2 //3 

makes sense as a bundle, with a well-defined bundle map ir[(t, e)] = [(£, 7r(e))] because of 
the equivariance. 

Recall how an automorphism acts on global sections: if s : S 1 x — > Ind2 is a global 
section, then 

~ -1 

s -> h o so h 

under the action of Is- Let us use the same notation for the lift of this action to the 
bundle of spinors. Now this action commutes with the Dirac operator, so we can ask 
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about the trace of 7 3 on the space of Dirac zero modes. In fact, we will be able to glean 
this information in a simpler way which also reveals the electrical charges of the states we 
have found. 

Note first that 7 3 defines a map 1^ which acts only on the 7 factor of Ind 2 in (3.12). 
Note that 7g squares to —1, the map which is —1 on each fiber but trivial on the base. 
In fact, prior to this section we have only dealt with 7. Global sections of Ind 2 have 
the form s(t,m) = (£, s(t,m)) G R x 7. The 5" 1 piece of the Dirac equation dictates the 
t— dependence of s to be s(t, m) = e tQt / 2 s(m), where s is a zero mode on 7 and Q is the 
electric charge. Sections of this form are clearly eigenstates of the electric charge operator 
— 2zjp Note that multiplication by a complex phase along the fibers is well-defined on 7 
since it commutes with transition functions and is therefore independent of trivialization 
(in a given trivialization, s takes the form of a complex- valued function). 

Now s will descend to a global section on Ind 2 if 

s(t : m) = I3 o s o 7 3 ~ 1 (t, m) 
= I 3 o s(t — 7r, I 3 ~ 1 m) 

= T 3 (t- 7T, e ^(*-)/2/W o s O 73-^) 

= {t,e- i ^' 2 e i ^ 2 li I) osoI- 1 m). 

We saw (I^) 2 = —1 on 7 and since 1^ acts on the space of zero modes we can take 
s to have definite eigenvalue equal to r (r = ±i). Putting 1^ o s o 7 3 ~ 1 m = rs(m) and 
reinstating the t— dependence yields 

r(-i) q = 1. (3.13) 

In fact, this analysis holds for any odd power of 7. For even powers, the bundle of zero 
modes is untwisted over the 5" 1 piece and therefore 1^ squares to the identity. Thus the 
eigenvalues r take the form r = ±1. The same condition (3.13) applies. Thus: 

n odd r = +i states have charge 1 mod 4 

r = -i " 3 mod 4 

n even =^r = +1 states have charge mod 4 

r = -1 " 2 mod 4. 

In other words, every zero mode of 7 yields a physical monopole solution with electric 
charge dependent upon its 7g eigenvalue. 
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How do we count the number of solutions with a given charge? We take the trace of 
I3 on the space of zero sections. This is done for n = 4 (the fourth power of I) below. 
For n = 3 we remark that there is only one eigenvector, with eigenvalue ±i, which must 
be the complex conjugate of the n = — 3 eigenvector. These states are charge conjugates 
(note that 1 and —1 have different odd values mod 4). This argument tells us nothing for 
the even states, but for n < 2 there are no solutions anyway. 

Let L(Is) denote the trace of 1^ restricted to the kernel of D on I 4 . In order to 
compute L(Is) we need a noncompact variant of the Atiyah-Segal-Singer equivariant index 
theorem. We may argue exactly as we did in section 3.2 to obtain the following proposition. 
(See m 



Proposition 3.6.1 

L(Is) =limL^oo(H'ms^o / trrI^e~ sD dx + / treox 

Jr<L JO Jr=L 

rl^De'^dads). 

The small s limit can be evaluated exactly as in the compact case. It localizes to a 
computation in a neighborhood of the fixed point set of I3, but ^3 has no fixed points. 
Hence (see, for example, Section 6.3 of fTF 



/imz / ^ 00 (/im s ^o / tTrljpe sD dx) = 0. 

Jr<L 

We are left to compute 

f°° f 2 
lirriL^oo / / treQTl^ De~ sD dads. 

Jo Jr=L 

We compute this defect term using the same semilocal parametrix we used in section 
3.5. The action of 1^ changes the computation in two ways. First, it acts by ±eie2 on 
the spinors. The sign is determined by the choice of the lift of the action of ^3 to the 
principal spin bundle. Secondly, 1% introduces introduces a rotation of it in the fiber of 
the fibration (3.9). This enters the computation by having us evaluate the parametrix not 
on the diagonal, but along the diagonal on the base of the fibration and on (ijj, i^ + tt) in the 
fiber. For the k th Fourier component, this is (— l) k times what one obtains by evaluating 
along the diagonal. So 1$ changes our computation by introducing fators of eie2(— l) fc . 
Once again trace identities make all terms vanish except 
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f rOO 

~ J2 tre 1 e 2 e 3 e 1 e 2 (-l) k e 3 (e L /2-l)(i[k + l]-(l/2-U 2 /4)e 1 e 2 )x 

JRP 2 £ JO 
e ^ Wfc+ l]-(l/2- f /V4)e 1 e 2 ) 2 (47rs) -3/2 ds 

/•oo 

= _ ( 27r ) V / tr(-l) fc (e L /2 - l)(iJfe - (1/2 - E^ei^e^"^-^/^^) 3 x 
(Ans)- 3 / 2 ds. 

This expression is skew under the involution that interchanges the eigenspaces of eie 2 and 
sends k —>■ —k. Hence the trace is zero. 

We can therefore conclude that for n = 4 

L(h) = 0, (3.14) 

and for n = —4 we also find L(I%) = 0. So there are an equal number of charge mod 4 
and 2 mod 4 two-monopole states. 



4. The BPS Spectrum 

Let us briefly review the predictions of Seiberg and Witten and compare our results 
with their predictions. Of particular interest to us are the global symmetry properties of 
the theory. It will be an important check that the BPS states we find appear in predicted 
representations of the global symmetry group. Let us first exhibit the flavor symmetry 
discussed in explicitly. For the purpose of determining the global symmetry, the La- 
grangian contains a kinetic term for the hypermultiplets given by the bilinear form 

M f M + M f M, 

where Nf flavor indices as well as the gauge group indices are suppressed. The chiral 
superfields M and M are in conjugate representations of the gauge group. However, 
since the gauge group is SU{2), the fundamental and anti-fundamental are isomorphic, 
so a symmetry mixing M and M is permitted. The other term to be preserved in the 
Lagrangian is the coupling to the Higgs field M T $M. If we define the 2iVj--dimensional 
complex vector 

V=( M + M-\ 
~ \i(M-M) J ' 
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then a symmetry V — > must preserve and 1/ T F i.e. A G U(2N f ) and A G 

0(2N f ; C). So At = A" 1 = A T , and thus A = A*; therefore A is in 0(2N f ) with real 
coefficients. At the quantum level, the parity in 0{2Nf) either reverses the sign of the 
electric charges (Nf odd), or is broken (Nf even). The relevant global symmetry group 
for states of a given charge is therefore SO(2Nf). 

As we discussed in section 2, the low-energy dynamics of monopoles and dyons is de- 
scribed by a supersymmetric sigma model. The Hilbert space decomposes into representa- 
tions of SO(2Nf). For magnetic charge k = 1, the bundle of zero modes is one-dimensional, 
and so the zero mode anti-commutation relations 

{ijj\^} = 6 ij i = l,...,2N f 

lead to spinorial representations of SO(2Nf) as noted in ||. For magnetic charge k = 2, 
the bundle is two-dimensional, and so the algebra that we must represent in terms of flavor 
properties is no longer a Clifford algebra but 

{^,^'t}=^, (4.1) 

where the ip 1 are now complex. This is just the usual annihilation and creation opera- 
tor algebra. BPS states with magnetic charge two therefore fall into representations of 
SO(2Nf) rather than its universal cover. 

For Nf = 0,1,2, we found no BPS bound states with magnetic charge two, and 
the singularity structure proposed by Seiberg and Witten required no bound states. In 
the case Nf = 0, the semi-classical BPS spectrum is completely determined since there 
are no bound states for any magnetic charge k > 1. Such states would correspond to 
normalizable anti-holomorphic forms on a non-compact Calabi-Yau manifold, and there 
are no such forms. 

For Nf = 3, we found a single bound state for each charge (2, n) where n is any odd 
integer. These states are singlets under the SO (6) flavor symmetry since they correspond 
to the Fock vacuum for the algebra (4.1) and its complex conjugate. In this case, the 
singularity structure of the moduli space can be analyzed by looking at the limit of three 
equal very massive quarks. Since the mass is large, there is a singularity in the semi- 
classical region of the moduli space where the scalar vacuum expectation value is large 
and the quarks (in the 3 of the SU(3) flavor symmetry) become massless. In the strong 
coupling region one can integrate out the massive quarks semiclassically and relate the 
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singularities to the two known singularities - a monopole with (n m ,n e ) = (1,0) and a 
dyon (1, 1) both in the 1 of SU (3) - of the Nf = quantum moduli space. In this way, the 
nonabelian global charges of the massless states at the singularities are calculated. These 
charges cannot change. One can then let the masses of the quarks go to zero and determine 
the representation theory under the full Spin(6) = SU(4) symmetry. This requires the 3 
and 1 to combine into a 4 singularity, while the other singularity (now a 1 under SU(4)) 
remains separate. For SU(4) the center Z4 is determined by general consistency conditions 
to act as e l7r ( n ™+ 2n e)/ 2 ? yielding the condition for the 4 that n m + 2n e = 1, with minimal 
solution (n m ,n e ) = (1,0). This just the standard (1,0) BPS monopole. The singlet state 
is trivial under the center and thus has minimal solution (n m , n e ) = (2, 1). As explained in 
0, this state should be continuously connected to a BPS state which exists semi-classically, 
and we have found such a state. 

The case of most interest is Nf = 4. This theory is conjectured by to be self- 
dual under SX(2,Z). However, the situation is somewhat different from the N=4 Yang- 
Mills theory which is also conjectured to be self-dual. For N=2, Nf = 4, the SL(2, Z) is 
believed to act on the representation spaces of the Spin(8) global symmetry as well. The 
elementary hypermultiplets transform in the vector of Spin(8) (i.e. the fundamental of 
SO (8)), while the (1,0) and (1, 1) monopole states transform in opposite chirality spinor 
representations. Although these representations are not isomorphic, they are all eight 
dimensional and permuted by the S3 which acts as outer automorphisms of Spin(8) (an 
inner automorphism would give an isomorphic representation). The trivial representation 
is trivial under this S3 as well. We will not further motivate this prediction, but will rather 
simply state the action of symmetry. The action of SL(2, Z) on a state is to transform 
the monopole numbers in the usual way - (n m ,n e ) — > (an m + bn e ,cn m + dn e ) under 



matrices mod 2 is easily seen to be isomorphic to S3) acting by left multiplication on the 
representations (0,0) = trivial o, (0,1) = vector v, (1,0) = spinor s, (1,1) = spinor c, 
where all numbers are defined mod 2. Simply stated: the representations are determined 
by (n m , n e ) mod 2. 

The BPS spectrum contains a stable elementary electron with charge (0, 1). The self- 
duality conjecture then implies the existence of bound states for all charges (p, q) where 
p and q are relatively prime integers. Specifically, bound states with charge (2, q) with q 




G SL(2, Z) - and the representation is transformed under p 6 S3 as above, where 
natrix y ^ J G SL(2, Z) with entries mod 2 (the group of unit determinant 
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odd must exist and appear in the vector representation of Spin(8). We indeed found such 
states from the bound state solution corresponding to the excitation of a single zero mode 
7*' | — 4 > on the vacuum with £7(1) charge —4, and its complex conjugate. The allowed 
electric charges are 4g + 1 for one state, and 4q + 3 for the other. 

The BPS spectrum also contains neutrally stable heavy gauge bosons of charge (0, 2). 
If such states exist as discrete states in the theory then we should expect to see their 
partners under SL(2, Z) with charges (2p, 2q) as discussed in 0. These states must all be 
singlets under Spin(8). Further, the heavy gauge bosons are part of BPS multiplets with 
spins < 1 unlike the electrons which are part of BPS multiplets with spins < |. From 
our computations, we have shown that four bosonic bound state solutions exist that are 
singlets under Spin(8). We found in (3.14) that the electric charges are 4g for two of the 
solutions and 4q + 2 for the remaining two. The existence of such states certainly implies 
that if the theory is self-dual then the heavy gauge bosons must exist as bound states at 
the threshold of decay into electrons. Two solutions are also required if one is to construct 
a BPS multiplet containing a vector particle. Our findings are certainly in accord with the 
proposed duality. 

To further support the supposition that the two solutions at n = ±4 with a given 
electric charge are members of the same BPS multiplet, we can examine the difference 
in fermion number between the different bound state solutions. We begin by noting that 
in the models under consideration, the fermion number is always integral. The fermion 
number of a bound state solution comes from two sources ||17|| : the first contribution is 
from fermions in the effective action for pure N=2 Yang-Mills. Since the moduli space is 
a product, the action can always be written as the sum of two terms: S e ff(k = 1) and 
an interacting piece Si n t- Quantizing the fermions from the first term gives us the usual 
four-dimensional BPS multiplet when acting on a spin zero vacuum. Our interest resides 
with the remaining fermions from Si n t, and the difference in the fermion number between 
the n = 3 and n = 4 bound state solutions. Viewing these fermions as spinors, and noting 
that the index has the same sign for n = 3 and n = 4, we can conclude that the bound 
state solutions at n = 3 and n = 4 have the same chirality. Therefore, the difference in 
fermion number from this source is zero mod 2. The other contribution, from the matter 
fermions, clearly produces a difference in fermion number. Therefore, we can conclude 
that the bound states at n = 4 differ in fermion number from the bound state at n = 3 
by one mod 2. Spin-statistics implies that the bound state at n = 3 is bosonic, and so the 
two states at n = 4 are fermionic. A BPS multiplet built on such a vacuum includes a 
vector particle as expected. 
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5. Conclusions 



We have studied the question of whether bound states of monopoles and dyons with 
magnetic charge two exist in supersymmetric Yang-Mills coupled to matter. This problem 
was solved by computing the number of L 2 solutions of the Dirac equation for bundle- valued 
spinors over the two-monopole moduli space. For Nf < 3, no bound states exist, while for 
Nf = 3, there is a single bound state for every odd value of the electric charge, which is 
a singlet under the SU(4) symmetry group. For Nf = 4, there is a bound state for each 
odd value of the electric charge. These bound states are in the vector representation of the 
Spin(8) flavor symmetry. There are also two bound states for each even value of the electric 
charge, which are singlets under the flavor group. Our findings provide dynamical evidence 
for the moduli space structure proposed by Seiberg and Witten - specifically, for the Nf = 4 
conjectured self-duality. To show the BPS spectrum of the Nf = 4 theory is truly SL(2, Z) 
invariant, similiar calculations are needed for higher magnetic charge. The main obstacle 
is the limited information about the metric for the higher charge monopole moduli spaces. 
Recently, the asymptotic metric on the k- monopole moduli space for the region where all k 



monopoles are far apart has been described ||18|1 . However, an understanding of the metric 
at the boundaries of codimension one would be desirable to extend computations of this 
type. 
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